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A CONCISE PARAMETRISATION OF AFFINE TRANSFORMATION 


SHIZUO KAJI AND HIROYUKI OCHIAI 

Abstract. Good parametrisations of affine transformations are essential to interpolation, de¬ 
formation, and analysis of shape, motion, and animation. It has been one of the central research 
topics in computer graphics. However, there is no single perfect method and each one has both 
advantages and disadvantages. In this paper, we propose a novel parametrisation of affine trans¬ 
formations, which is a generalisation to or an improvement of existing methods. Our method 
adds yet another choice to the existing toolbox and shows better performance in some appli¬ 
cations. A C-|—h implementation is available to make our framework ready to use in various 
applications. 


Throughout this paper all vectors should be considered as real column vectors, and hence, 
matrices act on them by the multiplication from the left. 

1. Introduction 

Affine transformation is an essential language for discussing shape and motion (see, for exam¬ 
ple, ID)- A common difficulty we often encounter while manipulating affine transformations is 
how to represent an element. A 3D affine transformation is represented by a 4-dimensional ho¬ 
mogeneous matrix (see eq. ([^), however, working directly with this representation is sometimes 
inconvenient since, for example; 

• the sum of two non-singular (or rotational) transformations is not always non-singular 
(or rotational) 

• various interpolation and optimisation techniques developed for Euclidean space do not 
apply straightforwardly. 

Many different parametrisations of transformations have been proposed which have certain 
good properties. Nevertheless, none of them is perfect and we have to choose one for each 
purpose. Here, we mean by parametrisation a map 0 from a space Vg to the set G of certain 
transformations. The idea is that one can operate on and analyse transformations in an easier 
space Vg instead of directly dealing with the complicated set G. 

The main goal of this paper is to propose a novel parametrisation of 3D affine transformations 
which possesses the following favourable properties: 

(I) Vg is a Euclidean space 

(II) 0 : Vg —!• G is continuous and differentiable 

(III) 0 is surjective and there exists a differentiable local inverse 0 : G —)■ Vg such that 
0(0(A)) = A for any A E G 

(IV) dim(VG) = dim(G) 

(V) For an important subclass of transformations H <Z G, there exists a linear subspace 
Vh C Vg with 0(1/^) = H 

(VI) 0 and 0 are computationally tractable at low cost 

The condition (I) makes available the basic operations on Euclidean spaces such as averag¬ 
ing as well as interpolation, calculus, and various linear analysis techniques such as Principal 
Component Analysis. Note that 0 converts linear operations to non-linear ones; addition in 
the parameter space Vg can be highly non-linear in G. The condition (II) is necessary if 
one wants to apply techniques from calculus such as differential equation to solve optimisation 
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problems including Inverse Kinematics. The condition (III) means that any transformation has 
a canonical representative in the parameter space and the correspondence is differentiable in 
both ways. The condition (IV) means that the transformation and the parameter has the same 
degree-of-freedom and there is no redundancy. The condition (V) means that the parametri- 
sation restricts to that of a subclass. For example, when we consider rigid transformations 
inside affine transformations, we can assure that an interpolation of rigid transformations is 
always rigid. The condition (VI) is mandatory, for example, for efficient creation of deformation 
animation, where millions times of computation of (j) and '0 per second are required. 

Note that there are many good parametrisations which do not satisfy one or more conditions 
listed above but come with other nice properties. We have to choose a parametrisation which 
suits a particular application, and the purpose of this paper is to add another choice to the 
existing toolbox. 

To clarify and demonstrate the meaning of the above conditions, we hrst look at the example 
of the 2D orthogonal group SO(2) consisting of two-dimensional rotations. For G = SO(2), we 
can take Vg to be the one-dimensional Euclidean space M, so that dimG = dimVc = 1. We 
take 0 to be the exponential map 

^ exp(«, ^ (-« —/) . 

where • Then 0 : Vd —)■ G gives a differentiable, surjective map, and it is also a 

group homomorphism, i.e., 4>{0i + 62 ) = 0 ( 0 i) 0 ( 02 )- The inverse map d • G —)■ Vd is given by 

if fa -b\\ _ ^ + kn ( 07 ^ 0 ) 

\\b a J J "^ 71/2 + kTT (a = 0,6 = ±l)’ 

where —ti/2 < tan“^ ^ < 7 r /2 is the principal value of the arctangent and k G Z. Note that 
d cannot be taken as a globally continuous map on G, but locally around any point of G it 
can be taken to be a differentiable map by choosing k E Z appropriately. Therefore, in this 
situation, the conditions (I), (II), (HI), and (IV) are satished. The condition (V) is out of 
question here since there is no interesting subclass of SO( 2 ). As for (VI), the computational 
cost of d and d is reasonable since we have the closed formulae as described above. When G is 
a more complicated class, we cannot expect this kind of simple solution. 

The mathematical tool we employ to construct our parametrisation is Lie theory, particularly, 
the Lie group-Lie algebra correspondence and the Cartan decomposition of Lie algebras (for 
mathematical background, we refer the reader to textbooks on Lie theory such as [M])- We 
also discuss several applications of our parametrisation. This paper is an improved and extended 
version of [H] . Introductory explanations of our method including the background can be found 
in 1111321 [33]. 


2. Related Work 

Precedent researches have given a number of different parametrisations to the subclasses of 
3D affine transformations listed in Figure]^ (see also (Q. One of the earliest example is the 
parametrisation of the group of 3D rotations SO(3) by the Euler angle. As is well-known, it 
suffers from the gimbal lock. In our language, the condition (III) is not satished. The gimbal 
lock issue is avoided when SO(3) is parametrised with the invertible quaternions. However, the 
dimension of the space of invertible quaternions is 4, which is greater by one than that of SO(3). 
This means we need an extra variable and the condition (IV) is not satished. The parameter 
space is almost Euclidean; except that the origin 0 of the parameter space does not correspond 
to any transformation. Although this is only a single point, it can be problematic. For example, 
we have to be careful when interpolate transformations so that the interpolation curve does not 
pass the point. If we restrict ourselves to the unit quaternions, we have the parameter space 
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Aff+(3) 

Sim’''(3) 

GL+(3) 


SE(3) 

CO+(3) 

Sym+(3) 

M3 

SO(3) 

M+ 


symbol 

transformation 

13 

the group of translations 

SO(3) 

the group of rotations (or the special orthogonal group) 

M+ 

the group of dilations 

SE(3) 

the group of rigid transformations (or the group of Euclidean motions or screw motions) 

CO+(3) 

the group of linear conformal transformations with positive determinants 

Sym+(3) 

the set of scale-shear transformations (or the set of positive dehnite symmetric matrices) 

Sim’''(3) 

the group of similarity transformations with positive determinants 

GL+(3) 

the group of linear transformations with positive determinants 

Aff+(3) 

the group of affine transformations with positive determinants 


Figure 1. Hierarchy of 3D transformations. The upper classes contain the lower classes. 


of the same dimension as SO(3). However, the parameter space is non-Euclidean and special 
techniques are required for basic operations (see, for example, [sn El EH El ESI El). 

The group of rigid transformations SE(3) can be parametrised with the dual quaternions, 
which is a generalisation of the parametrisation of SO(3) by the quaternions. In |2I]; they used 
the dual quaternions to blend rigid transformations and their method works particularly well 
in skinning. However, for other applications we may be troubled by the complicated structure 
of the parameter space; the space of the unit dual quaternions is the semi-direct product of the 
group of the unit quaternions and 

Let Aff^(3) be the group of 3D affine transformations with positive determinants. Note 
that an (invertible) affine transformation has the positive determinant if and only if it contains 
no flip (or reflection). In |1D] a method to interpolate elements of Aff’^(3) using the polar 
decomposition was introduced. Transformations are decomposed into the rotation, the scale- 
shear, and the translation parts, and then SLERP was used for interpolating the rotation part 
and the linear interpolation was used for the rest. This idea to look at Aff’''(3) as the (non- 
direct) product of three spaces SO(3), Sym’''(3), and has been fundamental and many of the 
current graphics systems adopt it. However, the parametrisations of SO(3) by the quaternions 
and Sym’'"(3) by matrices are not Euclidean. 

On the other hand, in |2] a dehnition of scalar multiple and addition in Aff^(3) is given 
based on the idea to parametrise Aff’''(3) by the corresponding Lie algebra, which gives a 
Euclidean parameter space. This is a generalisation of imi, where SO(3) is parametrised by 
its Lie algebra. The same idea is also used in [3H]. A notable feature of their construction, 
which is missing in jlO], is that the scalar multiplication satishes “associativity.” That is, for 
a, /d G M, the a-multiple of the /d-multiple of a transformation is equal to the a/3-multiple of 
it. However, a major defect of their construction is that it does not work with transformations 
with negative real eigenvalues. That is, there is no representatives for some transformations 
and the condition (HI) does not hold. This causes a big problem (see );4.4). Our work can be 
considered as a workaround of this inconvenience at the cost of loosing the associativity. In 
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addition, our parametrisation comes with a few advantages including a fast closed formula and 
better handling of large rotation. 


3. The parametrisation map 


The key idea of our method is to linearise the curved space of the Lie group of the affine 
transformations using the Lie algebra. This direction has already been pursued in the computer 
graphics community, for example, in uni El ES]. However, we have to be careful about the 
following points when we deal with the group Aff’''(3); 


the correspondence between the group Aff^(3) and its Lie algebra is not one-to-one. 
working in the Lie algebra is not intuitive since the meaning of each parameter is unclear 
(for example, finding consistent logarithm discussed in 14.3 becomes difficult), 
it involves the high cost computations of the matrix exponential and logarithm. 


To avoid those shortcomings, we propose a handy and mathematically rigorous parametrisa¬ 
tion. Using the Cartan decomposition of the Lie algebra of GL’'"(n), we establish a parametri¬ 
sation of Aff^(3) which satishes all the conditions listed in the previous section. (For (V), we 
can take all the important subclasses listed in the previous section as H.) 

Recall that Aff^(3) denotes the group of 3-dimensional affine transformations with positive 
determinants, i.e., the connected component including the identity. In other words, Aff^(3) is 
the group of 3D orientation preserving affine transformations. As usual, we represent elements 
in Aff"''(3) by 4 X 4-homogeneous matrices. 


( 1 ) 


Aff+(3) 


\ 

^Oii 

012 

Ol3 

l.\ 

] 

Li= 

021 

022 

023 

ly 

h 

det(A) > 0 / 


O 3 I 
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033 

1 

^0 

0 

0 

V 
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We often denote by the (i, j)-entry of A. From this representation, it is clear that Aff’''(3) 
is a 12-dimensional Lie group. We call the upper-left 3x3 part of A G Aff’''(3) as the linear 
part of A and denote it by A, and we denote by 'ipiiA) the translation part of A 


Let M{N) 
given by 


be the set of A^ x A^-matrices. 


/l 0 
0 1 
0 0 
\0 0 

Denote 


0 

0 ly 

1 h • 

0 ij 

by L the standard inclusion M(3) —)■ M(4) 


i(B) 




Then, A = '0i(A)t(A) for A G Aff’'“(3). 

We dehne our 12-dimensional parameter space as the product of two 3-dimensional and one 
6-dimensional vector spaces: 

^Aff+(3) := X 50(3) X sr)m(3), 

where 

so(3) := {X G M(3) | X = -W} 

is the set of the 3 x 3-anti symmetric matrices (which is the Lie algebra of SO(3)) and 

50m(3) := [Y G M(3) | T = V} 
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is the set of the 3 x 3-symmetric matrices. To sum up, our parameter space ^ 45 +( 3 ) is the 
12-dimensional Euclidean space of the form 


X so(3) X si)m(3) 


/ /l 0 0 li\ 

010/2 
001/3 
y \0 0 0 IJ 


0 

X 4 


/yi 

2/8 

2/9 

— X 4 

0 

X 6 

, Vs 

2/10 

2/11 

-X 5 

-xe 

0/ 

V//9 

2/11 

2/12 


\ 


where we identihed with the translation matrices so that the identity matrix serves as the 
origin. We emphasise that there is no restriction on the parameters and each variable can take 
any real number. 

The parametrisation map is dehned by 


(2) 0 : X so(3) X st)m(3) —)■ Aff’''(3) 

{L,X,Y) I—)■ L ■ i(exp(X) exp(F)), 


where exp is the matrix exponential. A locally differentiable inverse ip is dehned as follows: 

(3) 0 = : Aff+(3) -)■ X so(3) X sr)m(3) 

A ^ (0i(A),log(i5“^),log(^)), 


where S = \/^AA. That is, for any A G Aff’''(3) we have ip(A) G Vyjf+( 3 ) (unique up to modulo 
27r, as we see later in 14.3) such that 

0(0(A)) = A. 


It is obvious from the dehnition that our parametrisation satishes the conditions (I), (II), and 
(IV) in §1. For (V), we set 


^Sim+(3) — e ^Aff+(3) I C G M} 

W(3) = {a4,X,V)G VAff+(3)} 

^SE(3) = {(-^,-A, 0) G VAfj+( 3 )} 

^CO+(3) = {(-^4,-V, c/ 3 ) G VAfj+( 3 ) I C G M} 

^Sym+(3) = {(-^ 4 , 0 , V) G VAff+( 3 )} 

Mr 3 = {(-/>, 0,0) G VAff+( 3 )} 

^SO(3) = {(/ 4 ,-A, 0) G Vyfj+( 3 )} 

Vffi+ = {(/ 4 , 0, c/ 3 ) G Vyg+( 3 ) I c G M}, 

where J 3 ,14 are the identity matrices. When restricted to the above subspaces, our parametri- 
saion gives those for the corresponding subclasses of the transformations. We give explicit 
fomulae of 0 and 0 in the next section, and see that the conditions (III) and (VI) are satished 
as well. 


Remark 3.1. The product exp(X)exp(V) in ([^ is in fact the polar decomposition, where 
exp(X) G SO(3) is the rotation part and exp(V) G Sym’'“(3) is the scale-shear part. Therefore, 
there are a few choices for the order of the product in (|^. (Then, ([^ would change accordingly.) 
However, we believe the current order is the most intuitive when the user works directly in the 
parameter space. Usually, we are comfortable to think of a rigid transformation as a motion so 
that it is specihed with the global frame. On the other hand, scaling and shearing “changes” 
the shape of an object so it is more natural to think that it takes place with the local frame. 
Therefore, we would like to hrst perform scaling and shearing, which is parametrised by si)m(3), 
so that the standard frame of and the local frame are aligned. Then, the rigid transformation, 
parametrised by x so(3), is applied. 
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4. Computing the parametrisation map 


In this section, we give closed formulae for the parametrisation maps ([^ and ([^. For a C++ 
implementation, we refer the reader to |16j . 

The matrix exponential is defined by the infinite series as in Appendix]^ Since computation 
by the infinite series is very slow, it is crucial to have an efficient algorithm for applications. 
We give fast closed formulae including ones for the exponential and the logarithm of symmetric 
matrices, which have their own interests. 

4.1. Closed formula for ([^. First, we look at the definition (|^. The closed formulae for 
the exponential maps are given as follows. For 5 o( 3 ) part, Rodrigues’ formula [ 6 ] computes: 

( 4 ) 

cos^ 


exp(X) = Is + 


sin 9 

~r 


X + 


02 


-X" = h + sinc( 0 )X + - 


9 

smc- 

2 


X 2 


9 = 


tr(XX) 


When |0| is very small, the following second order approximation can be used to avoid small 
denominator: 


( 5 ) 


sinc( 0 ) 


i-C 

6 


For Y G si)m(3), the computation of exp(F) can be done by diagonalisation (see (0 in 
the appendix). However, we introduce a faster method avoiding the high-cost computation of 
diagonalisation (see Appendixj^for the derivation): Observe that the eigenvalues Ai, A 2 , A 3 are 
the roots of the characteristic polynomial and obtained by the cubic formula. This means, Aj’s 
can be computed by a closed formula and the computational cost is much cheaper than that of 
diagonalisation. Assume that Ai > A 2 > A 3 , and put Z = Y — A 2 / 3 . Then, the eigenvalues of 
Z are A^ = Ai — A 2 , X '2 = 0, Ag = A 3 — A 2 . Put 

'^1'^3(62(A'i) - e2(A3)) 


6 = 1 - 


1 

^~2 


A'l — Ag 

\'i{2e2{X'i) - 1) - A3(2e2(A3) - 1) 

2(a; - A') 


where 62 ( 0 ;) = 

( 6 ) 


exp(a:) — 1 — x 




. Then, we have 

exp(F) = exp(A 2)(/3 + bZ + cZ"^). 


We may have small denominators when some of Aj collide, and thus, A)^ —>■ 0 and/or Ag —)■ 0. 
Nevertheless, the above formula converges. In practice, when x is small we can use the second 
order approximation 

1 K x'l 

e 2 (A,)^- + - + 

4.2. Closed formula for (|^. To compute ([^, we provide formulae for log(S') and log(R), 
where S = \/^AA and R = 

First, note that *AA is symmetric positive definite, and hence, the square root and its loga¬ 
rithm are uniquely determined. We can calculate them by the diagonalisation as well, however, 
we adapt a similar method to the one given in the previous subsection. Denote the eigenvalues 
of ^AA by Ai, A 2 , Ag > 0. Assume that Ai > A 2 > Ag and set Z = MA/A 2 . The eigenvalues of Z 
are A/ = A 1 /A 2 ,1, Ag = A 3 /A 2 . Define 

A'£ 2 (A))-A;£ 2 (A' 3 ) 


CL — —1 -|- 


A'l- 

C2{X[)-C2{XA 

X[ - A'g 


A( 


c = 
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where £ 2 (t) = 

(7) 


log(a:) 


(X 


1 ) 


X 


Then, we have 


log(5') = ^ logfii) = i ((a + log(A 2))/3 - (a + c)Z + cZ^) 


.htl I 

O' Q • 


To obtain 


When A' —)■ 1, we use the second order approximation /^ 2 (A') - 
R, we compute S~^ = exp(—log(S')) by the formula ([^. Since the eigenvalues of — log(S') are 
— log(Ai)/2, we can reuse the values of A*. 

For R G SO(3), by inverting Rodrigues’ formula we have 

/Tr(R)-l' 


( 8 ) 


log(R) = 


-AR-^R), e = cos- 


{■ 


2sinc6' 

where cos“^ takes the principal value in [0,7r]. When 9 is small, the approximation ([^ can be 
used. When n — 9 is small, there is a different method. Let v = {vi,V 2 , v^) be the eigenvector of 

/ 0 -V 3 ^2 \ ^ /Tr(R) - 1 

R with eigenvalue one. Then log(i?) = 9 \ 0 —Vi , where 9 = e cos“^ I - 

V-Ts V, 0 J V 2 

{v2{R _* 77 ) 13 ^ 0 ) 

“ Alternatively, one can look at the diagonal entries in Q 


-*R)i3<0) 

[?, Appendix A. 1.1] for details) 


with e = . , , 

-1 {v2{R 

to compute log(i?) (see 

As we mentioned before, here we have indeterminacy of cos“^ up to modulo 27r. This reflects 

the fact that our parametrisation can handle (or distinguish) rotations by more than 27r. If we 

impose continuity, we will have one explicit choice as is described in the next subsection. 

4.3. Indeterminacy up to 27r. The dehnition of ([^ contains an ambiguity up to a factor of 
27r. In this section, we see how we give an explicit choice and how it is useful in some cases. 

Let R = G SO(3), which appears in the dehnition (3). One could simply take the 

principal value of the logarithm, i.e., choose 0 < 0 < tt in (^. On the other hand, it is 
sometimes convenient to distinguish the 27r-rotation around the x-axis from the identity, and 
also from 27r-rotation around the i/-axis. For example, if we consider a rotational motion rather 
than the hnal position of it, we have to track the rotational degree larger than 2ti. (see Figure]^ 
where the target shape is twisted over 27r degrees). 

In such a case, we can choose a particular value for the logarithm using the coherency or 
continuity regarding the metric on the parameter space. Algorithm computes the logarithm 
of R which is closest to a given X'. 


4.4. Computational Efficiency. We compare our closed formulae for exponential and loga¬ 
rithm given in with two widely used algorithms; Fade approximation and the scaling-and- 
squaring method by Higham f [T3l [14] ) implemented in Eigen library m version 3.2.8]), and 
diagonalisation ( [11] ). Our implementation is given in [16]. We computed the matrix exponen¬ 
tial and logarithm for randomly generated 10® matrices for 10 times and measured the average 
computational time. The timing is given in Table Error is measured by the maximum of the 
squared Frobenius norm \X — exp(log(X))||. among all the 10® matrices. 


Table 1. Timing: measured with 1.7Ghz Intel Core i7, 8GB memory, single thread 



exp of siim(3) 

log of Sym’''(3) 

Error in Sym'’“(3) 

Ours 

0.2402s 

0.2685s 

4.792 X 10“^^ 

Higham 

0.5454s 

5.930s 

6.356 X 10-29 

Diag 

0.5720s 

0.6195s 

3.419 X 10-23 


We also compare different parametrisation of Aff’''(3). Alexa/SAM directly computes log 
and exp for Aff’''(3) (suggested in [2l|36])- Polar-|-DQ computes the polar decomposition and 
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Algorithm 1: consistent logarithm of rotation matrices 


Input: rotation matrix R G S0(3), anti-symmetric matrix X' G 5o(3) 
Output; logarithm of R closest to X' 

begin 

e ^ cos-i((Tr(/2) - l)/2) 

if 0 = TT then 


if tr(^X'X')l2 > 0 then 
vr 


return 


else 


sJtr{X'X')/2 


X' 



2 sine 

-/2 


{R - R) 


e' ^ ^/x% + x'l, + x% if X12X' 

I X ^ -X, -9 


12 


+ Xi3X'i3+X23X'23 


end 

while 9' — 9 > TT do 
I 9^9 + 27 ^ 

end 

while 9 — 9' > TT do 

I 9 ^ 9 - 2 tt 

end 

return 9X 
end 


< 0 then 


then nse dnal qnaternions for the rotation and the translation parts (snggested in [21]), and 
positive dehnite symmetric matrices for the scale-shear part. Polar-l-Grassia also compntes the 
polar decomposition bnt uses log for the rotation part (suggested in jHj). For computing the 
polar decomposition, we used the algorithm in [?], which is more efficient than by SVD and 
widely used in the computer graphics community (see [10]). Error is measured by the maximum 
of the squared Frobenius norm |X — '^(0(X))||, among randomly generated 10® matrices with 
det(X) > 10“^. We noticed that for X with a very small determinant, the error tends to 
get bigger with our method. Hence, to deal with near singular matrices, it is better to use 
Higham’s algorithm to obtain the polar decomposition and then use our formulae for exp and 
log to compute the parametrisation. 


Table 2. Comparison of various parametrisations 



Surjectivity 

Extrapolation 

Associativity 

Large rotation 

b 

b 

Error 

Ours 

Yes 

non-singular 

No 

Yes 

0.4835s 

0.3543s 

1.936 X 10"^^ 

Alexa/SAM 

No 

can be singular 

Yes 

No 

12.88 

0.8432 

3.269 X 10^ 

Polar+DQ 

Yes 

can be singular 

No 

No 

0.5366 

0.03531 

7.025 X 10"^® 

Polar+Grassia 

Yes 

can be singular 

No 

Yes 

0.5777 

0.1053 

7.371 X 10"^® 


Surjectivity means that any element X G Aff’''(3) can be parametrised. Extrapolation indi¬ 
cates whether extrapolation of two elements with any weights can fall out of Aff’''(3) or not. 
Surjectivity and extrapolation affects robustness of the method. Associativity means the asso¬ 
ciativity of scalar multiple explained in ^ Large rotation indicates the capability of dealing 
with rotations greater than 27r (see ^4.3). 
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Figure 2. Upper row: The four strings (left) are being deformed by our method. 

Two probes are placed above and below the strings, one hxed to the identity and 
the other rotated and translated. Notice that more than 27r rotation is recognised. 

Lower row: The leg of the octopus model (left) is being deformed by a probe 
rotated by large degrees. With our method (middle) the leg is made swirl, while 
with quaternions (right) one sees almost no deformation. 

Notice that due to the lack of surjectivity, Alexa/SAM has a huge error (see ^for theoretical 
background). 


5. Application 

In this section, we will give a few applications of our parametrisation with comparison with 
existing techniques (see also mi for a further comparison). Source codes with the MIT licence 
are available at https://github.com/shizuo-kaji/. 

5.1. Blending transformations. One direct application of Euclidean parametrisation of trans¬ 
formation in general is blending different transformations (see, for example, [?]). Suppose that 
transformations {A* | 1 < ^ < n} and the corresponding weights {wi G M | 1 < z < n} are 
given. One way to blend them is to take the linear sum in the parameter space: 

/ n 

(9) Blend(wi,..., Wn, Ai,..., A„) := 0 I ^ Wii){Ai) 

\i=i 

Note that the above formula reduces to interpolation of two transformations when n = 2 and 
W 2 = 1 — Wi- One can also apply standard techniques such as B-spline for interpolating three 
or more transformations. 

With the condition (V), if A/s are all in the same class, the blended transformation (|^ 
always stays in the same class regardless of the weights. 



5.2. Shape deformer. As a simple application of the Blend function, we construct a shape 
deformer (see Figure]^ similar to the one developed in [26]. A much more elaborated version 
is discussed in IZD] using our parametrisation, so we just give a simple idea of it here. 

Assume that a set of affine transformations {A* G Aff’'’(3) | 1 < z < n} (specified by “probe” 
handles) is given. We want to deform a mesh with vertices U, according to a given weight 
function zuj : U — )■ M (1 < f < n). The weight can be painted manually or calculated 
automatically from the distance of the vertex and the probe. 

Using the blending function defined in §5.1 we compute the deformed position of zz G U by 

u HA- Blend(zni(zz),..., Wn{u), Ai,..., A„) • u. 


See [20| for the details. 
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5.3. Shape blender. The authors considered in [12] to apply our parametrisation to morph 
two isomorphic 2D meshes based on the idea in [3|. It was extended to a shape blending 
algorithm for an arbitrary number of 3D meshes in im. We briefly review a simple version of 
it. 

Suppose that we are given a rest mesh Uq and target meshes Ui (1 < i < n). We assume 
that all the meshes are compatibly triangulated, i.e., a one-to-one correspondence between 
triangles of Uq and each Ui is given. We want to produce a blended shape U{wi ,..., with 
respect to the user specihed weights {wi G M | 1 < i < n\. We require that U{wi,... ,Wn) 
interpolates the given shapes, more precisely, U{wi,... ,Wn) = f/o if tCj = 0 (1 < Vi < n), and 

U{wi,... ,Wn) = Uk if Wi = 

First, we associate for each face fij (1 < j < m) of Ui the unique affine transformation Aij 
which maps the corresponding face foj to fij and the unit normal vector of foj to that of fij. 
Then, we put 

...,Wn) := Blend(wi, ...,Wn, Aij ,..., Anj) G Aff+(3). 

Note that those blended transformations are not coherent on the edges so we cannot apply 
them directly to Uq. In order to obtain a blended shape, we have to “patch” A'-’s to obtain 
a set of affine transformations consistent on the edges of Uo, that is, a piecewise linear trans¬ 
formation. This is done by finding the minimiser of the error function • ■ ■ > — 

A'j{wi,... to obtain a piecewise linear transformation Aj{wi^... ,Wn) (1 < j < m). A 

blended shape is then computed by 

U{Wi,...,Wn) = IJ Aj{wi,...,Wn)foj 

(see im for the detail). This is useful, for instance, to produce variations of shapes from a 
given set of examples (see Figure [^. 

5.4. Pose interpolation. Our method also has a direct application to interpolating transfor¬ 
mations (cf. [221 ESI [36]). Given a set of key frames for positions and orientations of objects 
such as cameras or joints. That is, we are given for each object, a sequence of affine transforma¬ 
tions Aj (1 < j < m). We want to interpolate the key frames to produce an animated poses 
(see Figure]^. Since our parametrisation takes value in the 12-dimensional Euclidean space, 
any ordinary interpolation method such as Hermite polynomial, Bezier, and B-spline can be 
applied in the parameter space: Let Interpolate(t, {knots}) be an interpolation curve with a 
given set of knots in Then, we can compute interpolated poses by 

A{t) := 0(lnterpolate(t, {fj{Aj) | 1 < j < m})). 

With the condition (V), if Aj^s are all in the same class H, the interpolated transformation 
always stays in the same class H regardless of the value of the time parameter t. 

6. Discussions: Comparison to precursors 

In this section, we discuss the relationship between our method and previous techniques from 
a theoretical point of view. Existing parametrisations of transformation are classified roughly 
into three types: based on matrix, Clifford algebra, and Lie algebra. 

Using matrix, we can parametrise Aff''"(3) by 4 x 4-matrices as in ([^. It is easy to work with 
and computationally efficient. Also it is equipped with the additive structure (in addition to the 
multiplicative structure) so that transformations can be blended by taking the weighted sum. 
However, the main drawback is that there are matrices which do not correspond to elements in 
Aff"''(3), namely those with vanishing determinant. This is problematic in many applications. 
For example, addition and scalar multiplication are not closed in Aff’'"(3). If you simply take the 
weighted sum of matrices, you may get degenerate transformation. Furthermore, the operations 


1 {i = k) 
0 {i ^ k) 
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Upper row left to right: rest shape Uq and target 
meshes Ui,U 2 - Lower row: shapes are blended with 
Wi = W 2 = 0.8 by onr method (left), by [21 ES] (mid¬ 
dle), and by polar decomposition and qnaternions [iU] 
(right). The left pictnre looks the most natural. The 
thorn in the middle picture is slanted because [21 EH] 
mixes up rotation and shear parts. In the right picture, 
the thorn totally collapse because quaternions cannot 
handle large rotation. 


ir 



Three yellow penguin shapes are 
blended with our method to produce 
variations (white). 


Figure 3. Shape blending 




Figure 4. Pose Interpolation: red rectangles are interpolated by the B-spline 
curve in the parameter space. Left to Right: our method, matrix logarithm [2], 
and homogeneous matrices (linear). The linear method shows degeneracy. 


do not restrict to the important subclasses. For example, the sum of two rotational matrices is 
not rotational. 

Quaternions (respectively, dual quaternions ([21])) can be used to parametrise SO(3) (re¬ 
spectively, SE(3)). More generally, using Clifford algebras one can parametrise various kinds of 
transformations. For example, the anti-commutative dual complex numbers ([2H]) parametrise 
SE(2) and conformal geometric algebra (CGA, for short) can be used to present Sim’^(3). In 
fact, CGA deals with a larger class of transformation including non-linear ones ([ElEDlSnillHl 
miEiss]). Glifford algebra based parametrisations are usually very efficient and the weighted 
sum can be used to blend transformations. However, there are two main disadvantages: hrst, 
it requires special construction for each class of transformation and, in particular, there are 
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no known construction for parametrising the entire Aff^(3). Among those listed in Figure]^ 
constructions for SO(3), SE(3), and Sim^(3) are known. The other disadvantage is that the 
parameter space is not Euclidean and one cannot directly apply techniques from linear analysis 
and calculus. In addition, the parameter space generally has greater dimension than the group 
of transformations to be parametrised. 

To cover the entire Aff’'"(3), Shoemake ([10]) suggested to hrst decompose an element of 
Aff'''(3) into the rotation, the scale-shear, and the translation parts by the polar decomposition 
and then parametrise rotational part by quaternions and the rest by matrices. This is one of the 
most commonly used method. However, the parameter space is not Euclidean; the rotational 
part is parametrised by the unit quaternions and the scale-shear part by the positive dehnite 
symmetric matrices. None of these parameter spaces is Euclidean and it causes problems in 
some applications (see Figure]^ and |^. 

In mathematics and physics, it is well-known that Lie algebras are useful to parametrise Lie 
groups. The exponential map is a differentiable map from the Lie algebra 0 to the Lie group 
G. When G is compact, this is surjective and gives a parametrisation of G by a Euclidean 
space 0 of the same dimension as G, which satishes (II) and (IV) discussed in For a Lie 
subgroup H of G, there exists a Lie sub-algebra f) and the image of the exponential map 
is contained in H. Hence, (V) is also satished. Grassia ([IH]) introduced the machinery to 
computer graphics community by suggesting to use the Lie algebra so (3) to parametrise SO(3). 
In [H], they combined this idea with the polar decomposition to blend elements in GL’''(3). 
They parametrised the rotation part by so(3) and the scale-shear part by Sym’''(3). 

In PEE], they used the Lie algebra to parametrise the whole Aff’''(3). However, Aff’''(3) is 
not compact and the exponential map is not surjective ((HI) fails to be satished). There exist 
transformations which are not parametrised; G G GL’'"(n) is in the image of the exponential 
map if and only ii G = for some B. In other words, log(A) for A G Aff''"(3) exists only 
when A = B'^ for some B (see [ 8 ], for detail). Although those transformations are not so many 
as was discussed in [36], in practice, it is problematic if one has to always take care of those 
exceptions (see Table 4.4[ ) and it does affect the result in some applications (see Figure [^. 
Another disadvantage is the cost of computation since the computation of the exponential and 
the logarithm usually involves iteration (see Table 4.4). 

Our construction is also based on the Lie algebra. To remedy the problem mentioned above 
we uses the Gartan decomposition in the Lie algebra, which corresponds to the polar decom¬ 
position in the Lie group. More precisely, for A G Aff^(3), exp (■^/j (A)) exp (-05 (A)) is the 
polar decomposition of A. This enables us to achieve (HI) while keeping the advantages of 
Lie algebra based methods. Note that thanks to our closed formulae for the exponential and 
the logarithm needed for our parametrisation, we avoid the cost inefficient computation of the 
polar decomposition. In fact, the above relation provides an alternative method to compute 
the polar decomposition. 

There are two main drawbacks of our method. First, with our parametrisation, the Blend 
function dehned in 15.1 fails to be bi-invariant, that is, invariant under the left and the right 
translation, whereas DLB developed in [21] for SE(3) has this property (see also [31], for 
interpolation). Bi-invariance is very important in some applications such as skinning. In [?], 
a bi-invariant mean of elements in a connected Lie group in general is studied. However, 
computing bi-invariant means and interpolation is very costly, and as far as the authors are 
aware, there is no parametrisation for Aff^(3) which achieves bi-invariance with a simple blend 
function. (In [2], they mistakenly claim that their method achieves bi-invariance.) Second, our 
Blend function lacks in associativity, which is satished by [21136]. 


7. Conclusions and future work 

We proposed a concise 12-dimensional Euclidean parametrisation for three-dimensional affine 
transformations. Our parametrisation has some good properties compared to the existing ones; 
most notably, it parametrises the entire Aff’''(3) by the Euclidean space of the same dimension 
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and it handles transformations containing large rotation. The Enclidean natnre allows us to 
apply various techniques developed for the Euclidean space. The computational cost is relatively 
cheap and a C++ implementation is provided. On the other hand, scalar multiple and addition 
in the parameter space do not map nicely to the transformation space; a straight line in the 
parameter space does not map to a geodesic in the transformation space. 

Since standard calculus and linear algebra are valid on our parameter space, we can apply 
differential equations, linear data analysis, and optimisation techniques on the group of affine 
transformations. Possible applications include inverse kinematics (|1I]), hltering captured data 
([25]), motion analysis and synthesis (|12]), and interpolation (|3l])- We will discuss them 
elsewhere. 

It would also be interesting to extend our result to parametrise wider classes of transforma¬ 
tions such as the projective transformations. 


Appendix A. Exponential and logarithm maps 

We recall the dehnition of the matrix exponential and the logarithm maps. For a square 
matrix A, the exponential map is dehned by 

(10) exp(/l) = f;T. 

1=0 

The “inverse” of the exponential is called the logarithm. It should satisfy the equation exp(log(A)) 
X. Let / denote exp or log. We have f{PAP~^) = Pf{A)P~^ for any invertible ma¬ 
trix P. Hence, when A is diagonalisable, that is, there is an invertible matrix P such that 
P~^AP = diag((ii, ^ 2 , • • •, dn), we can compute 

/(^) = /(^diag(di, (i2, ■ ■ ■, dn)P~^) 

( 11 ) = Pf{(Amg{di,d2,... ,dn))P~^ 

= Fdiag(/(di),/(d 2 ),...,/(dn))P-'. 

Note that log(A) does not always exist and even when it exists, it is not necessarily unique in 
general ([ 8 |). 


Appendix B. Derivation of the exponential and logarithm formulae in l|4] 


For Y G 5t)m.(3) we consider /(E), where / = exp. (And respectively, for Y G Sym’'“(3), we 
consider f{Y), where / = log.) The main idea is to divide the inhnite series (10) by the degree 
three characteristic polynomial to reduce it to a degree two polynomial (see, for example, |29j). 
More explicitly, let Py{x) = det(a ;/3 — Y) be the characteristic polynomial of Y. Then by 
dividing by the degree three polynomial py{x), we have 


f{x) = q{x)pY{x) +r{x), 

where r{x) = a + bx-\-cx'^ is the reminder. By Cayley-Hamilton’s theorem, we know Py{Y) = 0, 
so 


(12) f(Y) = r(Y) = ah + hY + cY\ 

Since Y is symmetric, we can diagnalise it to have Y = PDP~^, where D = diag(Ai, A 2 , A 3 ). 
To determine a,b,cE M, consider f{D) = diag(/(Ai),/(A 2 ),/(A 3 )) and 

f{D) = p-^f{Y)P = P-\ah + bY + cY^)P 

= ah + bD + cD^. 


Then, a, 6 , and c are obtained by solving the following Vandermonde’s linear system 


/I Ai An 
1 A 2 An 
yi A 3 \i) 
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Proposition B.l. When the eigenvalues are pairwise distinct, the coefficients in (12) is given 
by 

■s = /(Ai)/((Ai — A2 )(Ai — As)), 

t = /(A2)/((A2-A3)(A2-Ai)), 

^ = /(A3)/((A3 — Ai)(A3 — A2)), 
a = SA2A3 + tX^Xi + UX1X2, 
h = —s(A 2 + A3) — ^(As + Ai) — m(Ai + A2), 
c = s + t + u. 


(13) 


Next, we make this formula more robust so that it is valid even when some of the eigenvalues 
coincide. Assume Ai > A 2 > A 3 and put Y = Y — X 2 I 3 and f{x) = f{x + A 2 ). By this 
substitution, the computation of f{Y) is reduced to that of f{Y), where the eigenvalues of Y 
are Ai = Ai — A 2 , A 2 = A 2 — A 2 = 0, A 3 = A 3 — A 2 . Denote by T 2 {x) the analytic function 
f{x) - /(O) - f'{0)x 


x^ 


s 

t 

it 

a 

b 


By (13), the coefficients in f(Y) = al^ + bY + cY"^ are computed as 

/(Ai)/(Ai(Ai-A3)), 

/(0)/(A3Ai), 

/(A3)/((A3 — Ai)A3), 

tA 3 Ai = /(o), 

—sAs — t(X3 + Ai) — uXi 
1 

AiA 3 (Ai — A 3 ) 

1 


{“/(A^Ag — /(0)(Ai — A 3 )(A 3 + Ai) + f{X3)X^} 


{-(/(A0-/(0))A^ + (/(A3)-/(0))Af} 


AiA 3 (Ai — A 3 ) 

AiA3(T2(Ai)-T2(A3)) 


= f(0)- 


c = 


s + t + u 
1 


AiA3(Ai 

1 


A3) 


Ai - A3 

{/(Ai)A 3 + /(0)(Ai-A3)-/(A3)Ai} 
{(/(A0-/(0))A3-(/(A3)-/(0))Ai} 


AiA3(Ai — A3) 

AiT2(Ai)-A3T2(A3) 

Ai — A 3 

Although b and c may suffer from small denominators when Ai —)■ A3, it can be stably computed 
by substituting T 2 (Ai) with its Taylor expansion T 2 (Ai) = T 2 (A 3 ) + (Ai — A 3 )T 3 (Ai — A 3 ), where 
T 3 is an analytic function. Notice that since we assumed Ai > A 2 > A 3 , |Ai — A 3 I —)■ 0 implies 
Ai, A 3 —)■ 0. For / = exp we have 


exp(a;)-l-a; n 

T2 {x) = exp(A 2 )--= exp(A 2 ) 


x^ 


1 X X' 

2^3! 


to obtain 


b = exp(A 2 ) ( 1 
c = exp(A 2 ) 


A 1 A 3 

llT 


+ 


1 A 1 + A 3 A2 + AiA3 + Ai 

2^3!^ 4! 


These sequences converge quickly when Ai, A 3 —)• 0. 
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To derive the formula for log in 
log(F) - log(A2)/3. 


U.2 


we put Y 


Y /\2 and argue similarly using log(F) = 
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